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1 .0  \ 1 Introduction 


The  purpose  of  this  report  is  to  present  the  climatology  of  subsurface 
temperature  in  the  central  Pacific  Ocean.  This  climatology  provides  an  in-depth. 

'fine  scale*  description  of  the  thermal  structure  as  opposed  to  the  general  climatology 
developed  by  Robinson  (1976)  in  that  it  (a)  presents  finer  horizontal  and  vertical 
scale  information  on  the  temperature  field,  (b)  utilizes  data  not  included  in  previous 
climatologies,  and  (c)  represents  actual  observations  as  opposed  to  smoother  fields 
derived  from  actual  observations.  These  improvements  in  the  climatology  will  be 
useful  in  the  NORPAX  project  for  determining  the  depth  to  which  temperature  anomalies 
penetrate  in  the  ocean.  Also,  the  data  from  which  the  climatology  is  developed  can 
be  used  to  study  the  interannual  changes  in  heat  content  of  the  central  ocean.  It 
was  with  these  ideas  in  mind  that  the  present  work  was  undertaken. 


2.0  Data  Sources/Densitv 


The  initial  data  were  supplied  by  the  NORPAX  Data  Group  on  7-track,  800  BPS, 
character  magnetic  tape,  card  image  records,  and  multiple  records  per  BT.  There 
were  5 XBT,  19  MBT,  and  2 hydrographic  data  tapes  that  represented  all  of  the 
temperature/depth  in  a region  bounded  by  latitudes  25°N  to  60°N  and  longitudes  140°W 
to  180°W  for  the  time  period  1942-1974. 

Each  input  record  from  the  raw  data  tapes  was  tested  for  validity  and  data 
range  on  the  IBM-1800  computer  and  then  merged  into  a standard  header/data  format 
on  a group  of  4 working  tapes.  The  data  available  taken  from  each  observation  were: 


I type 


- Lat-Lon  Index:  latitude  = (IBX=3),  longitude  = (IBX  -1)*4 

- 1 = XBT,  2 = MBT,  3 = Hydro  for  observing  system 

- 25°N  to  (but  not  including)  60“N  in  degrees 

- Fractional  latitude,  minutes 

- 140°W  to  (but  not  including)  180®W  in  degrees 
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Lonin 

Year 

Month 

Day 

Hour 

Minute 

Reft 

Temperature 


- Fractional  longitude,  minutes 
= 1942  - 1974 

-1-12 

-1-31 

-1-24 

-1-60 

- Calibration  temp  (negative  for  no  data) 

- Data  0 0,  30,  60,  90,  120,  150,  200,  250,  300,  400, 
500  meters  (see  below) 


If,  In  processing  the  data,  a surface  temperature  was  missing,  and  there  was  a 
temperature  at  10  meters  or  less,  this  temperature  was  used  as  a surface  temperature. 

Linear  Interpolation  was  used  to  calculate  temperatures  at  standard  depths. 

The  Individual  types  of  Instruments  used  to  measure  the  temperature  profile 
are  discussed  briefly  below  (see  Wert,  1976  for  additional  Information). 

1 

Expendable  BT  (xBT) 

Data  are  sent  to  FNWC  In  the  form  of  strip  charts  and  their  accompanying  log  L' 

sheets.  These  data  are  usually  processed  on  the  semi-automatic  Calma  (model  403) 
digitizer.  The  Calma  device  records  the  minimum  number  of  points  necessary  to  re- 
construct the  original  trace.  Sometimes  an  analog  digital  data  system  Is  also  used. 

Regardless  of  the  digitizing  device,  the  data  are  recorded  on  magnetic  tape.  The  . 

repeatability  of  the  Calma  digitization  Is  stated  to  be  better  than  96X.  There  were 
12,639  XBTs  available  for  analysis. 

Mechanical  BT  (MBT) 

Data  are  from  the  National  Oceanographic  Data  Center.  The  data  base  Is  made 
up  of  digitized  records  obtained  from  MBT  slides  and  their  appropriate  log  sheets. 

Most  of  the  MBT  digitization  was  performed  on  contract  by  Scripps  Institution  of 
Oceanography  with  a semi-automatic  digitizer.  In  this  system,  both  the  grid  and  trace. 
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aligned  according  to  the  temperature  and  depth  correction,  are  projected  onto  a 
screen.  The  trace  Is  followed  manually  with  a stylus  and  the  temperature  data 
suitably  digitized.  Ancillary  Information  Is  simultaneously  recorded  directly  on  “i 

magnetic  tape  along  with  the  digitized  trace  data.  This  digitizing  method  has  the 
advantages  of  being  faster  and  somewhat  more  accurate  than  reading  by  eye.  In  all, 
there  were  57,328  MBTs  available  for  analysis. 

Hydrocast  Data 

All  of  the  hydrocast  data  are  archived  by  NODC  and  available  on  magnetic  tape. 

The  data  used  In  this  report  comes  from  the  1974  update  of  the  NODC  file  and  provided 
7,422  hydrocasts. 

Combining  all  data  sources  gave  77,389  (raw)  observations  of  temperature  vs 
depth  In  the  region  described  above.  The  distribution  of  these  data  throughout  the 
historical  record  (Figure  1)  shows  the  number  of  observations  per  month  In  the  entire 
study  region.  This  distribution  Is  rather  uniform  In  time  although  a higher  observe- 
tion  density  exists  during  summer  months.  Also,  note  the  phaseout  of  the  MBT  In  j 

1966  - 68  and  the  subsequent  phase-in  of  the  XBT.  Since  these  Instruments  measure 
to  different  depths,  one  may  expect  the  climatology  for  the  greater  depth  (>150m)  to 
be  representative  of  a shorter  base  period.  This  may  be  important  for  some  uses. 

Finally,  It  should  be  noted  that  ocean  weather  ships  contribute,  on  the  average, 
slightly  over  100  observations  per  month  to  the  total  data  set. 


3.0  Data  Editing 

Review  of  the  space/time  density  of  available  data  suggested  that  an  appropriate 
grid  for  construction  of  the  mean  field  would  be  quadrangles  2°  of  latitude  x 10®  of 
longitude.  The  time  Interval  for  the  climatology  was  selected  to  be  1 month.  The  ; 
latter  value  was  small  compared  with  the  time  scale  of  change  of  the  anomalies  which 
It  was  desired  to  study.  Similarly,  the  spatial  dimensions  of  the  averaging  areas 
were  also  small  compared  to  the  spatial  scale  of  anomalies.  Perhaps  as  Important, 
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the  small  latitudinal  extent  of  the  averaging  areas  (called  'boxes')  allowed  us  to 
avoid  the  spatial  aliasing  problem  that  results  from  taking  quasi  randomly  distributed 
samples  In  a high  gradient  field.  Thus,  the  fact  that  the  sampling  centroids  for  a 
given  box/month  do  not  occur  precisely  In  the  middle  of  the  box,  will  Introduce 
negligible  error  to  estimation  of  the  mean  field. 

After  all  of  the  available  BT  and  hydrocast  Information  had  been  sorted  Into 
appropriate  geographical  boxes  and  one-month  intervals,  it  was  gridded  In  the  vertical. 
The  standard  depths  used  were  0,  30,  60,  90,  120,  150,  200,  250,  300,  400  and  500 
meters.  This  selection  of  depths  allows,  us  to  adequately  resolve  the  seasonal  change 
In  the  heat  content  of  selected  slabs  of  ocean.  Note  that  it  would  be  somewhat 
difficult  with  the  coarse  grid  in  the  near  surface  region  to  determine  the  "depth  of 
the  mixed  layer,"  a nebulous  quantity  at  best. 

The  data  field  resulting  from  the  above  operations  was  a series  of  gridded 
values  in  x,  y,  z and  t space.  These  values  were  next  subjected  to  two  screening 
procedures  to  remove  data  that  was  clearly  incompatible  with  the  main  body  of 
information. 

Screening  Procedure  #1 

Climatological  values  of  temperature  at  0 and  120  meters  for  each  grid 
location  and  month  were  interpolated  from  the  Robinson  atlas.  These  values,  which 
had  been  carefully  scrutinized  and  smoothed,  presented  a first  guess  at  the  mean  field. 
The  individual,  raw  BT  information  was  compared  with  this  first  guess  at  the  two 
depths  indicated  above.  An  individual  BT  was  rejected  if  it  differed  from  the  first 
guess  by  more  than  ± 40^,  a value  representing  very  large  interannual  variations 
(Ballls,  1973a,  b,  c).  Thus,  any  BT  trace  that  failed  this  first  screen  was  eliminated 
from  the  data  set  altogether. 

Screening  Procedure  #2 

The  remaining  data  were  used  to  form  monthly  means  on  the  x,  y,  z grid  system. 
The  standard  deviations  about  these  means  were  also  computed.  Each  observation 
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belonging  to  a given  grid  point  and  time  was  checked  to  see  that  it  was  within  three 
standard  deviations  of  the  newly  computed  mean.  Values  outside  of  three  standard 
deviations  again  caused  the  entire  BT  to  be  excluded  from  the  data  set.  The  second 
screening  procedure  eliminated  only  a small  fraction  of  the  observations.  However, 
once  these  had  been  eliminated  a new  mean  and  standard  deviation  were  computed. 

The  means,  standard  deviation  and  number  of  observations  per  space/time  grid 
point  remaining  after  the  original  data  had  passed  through  the  two  screening  procedures 
represent  the  substance  of  this  atlas.  In  all  there  are  71,295  individual  observa- 
tions of  temperature  versus  depth  that  have  gone  into  the  establishment  of  mean 
fields  shown  in  Section  5.  Approximately  40%  of  these  observations  come  from 
standard  weathership  locations  and  so  the  quality  of  the  mean  field  estimates  is 
quite  spatially  dependent.  Nevertheless,  the  total  field  turned  out  to  be  rather 
smooth  and  very  'reasonable'  (see  below). 


4.0  Verification 

A number  of  methods  have  been  used  to  verify  that  the  resulting  estimate  of 
the  climatological  thermal  structure  field  in  the  central  Pacific  are  reasonable. 
These  are  enumerated  briefly  below. 

(A)  Vertical  consistency 

Over  much  of  the  region  of  interest  one  would  expect  the  temperature  to 
decrease  as  a function  of  depth.  The  exception  to  this  lies  in  the  region  near 
latitude  45°N  where  there  is  a well  known  thermal  inversion.  With  the  exception  of 
this  region,  all  other  data  were  inspected  to  insure  that  temperature  did  decrease 
with  increasing  depth. 

(B)  Horizontal  consistency 

One  would  expect  the  temperature  to  increase  as  latitude  decreases.  This 
was  numerically  checked.  With  the  vast  majority  of  the  data  this  is  true.  However, 
there  are  examples,  particularly  in  the  northern  part  of  the  study  area,  where 
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temperatures  at  higher  latitudes  during  a selected  month  are  slightly  higher  than 
those  Immediately  to  the  south.  This  situation  occurs  where  there  are  a very  small 
number  of  samples  available  to  estimate  the  means  In  a given  space/time  box.  These 
situations  have  not  been  excluded  from  the  tables  shown  in  Section  5.  Thus,  they 
can  provide  the  critical  reader  with  an  Idea  of  the  uncertainty  that  may  be  Inherent 
in  our  climatological  field. 

(C)  Temporal  consistency 

The  annual  cycle  of  the  surface  and  90  m mean  field  in  each  of  the  given 
spatial  boxes  has  been  plotted  and  inspected  to  insure  that  values  obtained  during, 
say,  summer  are,  in  fact,  greater  than  those  obtained  during  winter.  Examples  of 
these  plots  for  selected  latitude  bands  and  depths  are  given  in  Figure  2. 

(D)  Comparison  with  other  data  sources 

Ideally  we  would  like  to  compare  our  observations  with  the  atlas  results 
of  Robinson  (1976).  However,  that  data  source  has  already  been  used  in  our  first 
screening  procedure  and,  therefore,  such  a check  would  not  be  independent.  It  is 
possible,  however,  to  compare  the  surface  climatology  with  the  estimates  of  the  mean 
sea-surface  temperature  field  obtained  from  Jerome  Namias.  The  comparison  is 
presented  in  two  ways  in  Figures  3 and  4.  The  first  case  shows  a scatter  plot  of 
SST  determined  from  millions  of  ship  injection  temperatures  (Namias)  versus  those 
obtained  from  the  surface  values  of  the  more  limited  BT/hydro  data  set.  The  ship 
Injection  temperatures  are  seen  to  be  warmer,  on  the  average,  than  those  obtained 
from  the  BTs,  but  this  effect  was  expected  and  has  been  previously  documented  (e.g., 
Robinson,  1976;  Saur,  1963).  The  scatter  diagram  does  have  an 
interesting  property  in  that  the  two  estimates  of  SST  converge  where  the  highest 
temperature  is  observed.  Comparison  with  the  ship  injection  temperatures  may  be 
presented  In  another  way.  Figure  4 shows  the  frequency  distribution  of  the 
difference  between  the  two  SST  fields;  one  obtained  from  the  BTs,  the  other  obtained 
from  the  ship  Injection  temperatures.  It  Is  seen  that  the  mean  difference  between 
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them  is  approximately  0.4°C,  as  expected.  Taking  this  into  account,  it  is  abundantly  I 

clear  that  the  average  difference  between  the  two  fields  is  typically  less  than  a | 

degree  centigrade  and  that  the  difference  is  approximately  normally  distributed.  On 
the  basis  of  the  two  above  curves  we  may  conclude  that  the  SST  field  derived  from  ■ 

the  BT  information  is  statistically  equivalent  to  that  derived  from  ship  injection  | 

temperatures  (after  allowance  is  made  for  the  fact  that  the  injection  temperatures  j 

are  known  to  be  higher  than  the  bucket  temperatures  associated  with  the  BT  observations).  \ 
(E)  At  least  three  different  types  of  instruments  have  been  used  to  obtain 
the  data  that  has  gone  into  the  mean  field.  One  might  well  ask,  do  these  instruments 
bias  the  quality  of  measurement;  that  is,  do  the  mean  fields  obtainable  from  the 
different  Instruments  agree  in  a statistical  sense?  This  is  essentially  the  Berhens- 
Fisher  problem  and  can  be  answered  with  a version  of  the  Students'  't-test.'  Un- 
fortunately the  problem  is  compounded  here  by  the  fact  that  the  time  interval  covered 
by  the  mechanical  BTs  is  essentially  not  the  same  time  interval  sampled  by  the 
majority  of  the  XBTs  (cf  Figure  1).  Thus  the  instruments  could  be  returning  the 
same  estimates  of  a field  that  does  not  have  a stationary  mean  value  since  the  period 
of  the  '50s  and  early  '60s  may  not  be  similar  to  the  period  of  the  '60s  and  early  '70s. 
Hence,  the  t-test  mentioned  above  not  only  compares  the  quality  of  measurement  from 
the  various  instrument  types,  it  also  says  something  about  the  long  term  stability  of 
the  climatological  mean  field  developed  in  this  work.  There  is  a further  restriction 
that  the  mechanical  BTs  seldom  go  below  150  meters.  Hence,  any  comparisons  made 
between  the  different  types  of  instruments  must  be  made  at  this  depth  or  shallower. 

All  combinations  of  pairs  of  fields  have  been  examined  with  the  test  indicated 
above.  Also,  mean  field  derived  from  each  measurement  device  has  been  compared  with 
the  grand  mean  field  developed  from  all  of  the  data.  The  general  conclusion  derived  ; 

from  these  numerous  analyses  is  that  the  instruments  do  not  introduce  a significant  | 

bias  in  estimates  of  the  mean  field,  except  perhaps  at  the  sea  surface  where  the  | 

i 

comparisons  are  marginal.  1 
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5.0  Results/Cautions  | 

'j 

Tables  1 through  11  show  the  estimated  mean  field  for  each  month  for  each  j 

2 X 10  degree  quadrangle  and  each  depth.  Selecting  a given  time  and  depth,  each  box  j 

contains  three  entries.  The  center  value  gives  the  mean  temperature  In  degrees  j 

1 

centigrade,  the  value  on  the  lower  right  the  standard  deviation  associated  with  that  | 

mean  and  the  number  on  the  lower  left  In  the  box  gives  the  number  of  observations  J 

available  from  which  the  mean  and  standard  deviation  were  constructed.  The  latitude/  - ] 

longitude  coordinates  of  the  data  field  are  given  along  the  sides  of  the  printout  and  ] 

'i 

Indicate  the  centroid  of  the  sampling  regions.  All  of  the  data  contained  In  the  fields 
Is  available  on  magnetic  tape. 

ii 

The  resulting  mean  fields  do  have  several  peculiarities  that  require  caution  j 

on  the  part  of  the  user.  The  most  notable  events  are:  j 

(A)  In  some  cases  water  temperatures  are  higher  at  more  northerly  latitudes  than  , 

they  are  at  latitudes  Immediately  to  the  south.  This  event  occurs  very  Infrequently  and,  | 

when  it  does,  is  generally  only  apparent  at  the  highest  latitudes.  It  can  be  traced  | 

i 

directly  to  the  fact  that  there  Is  very  limited  data  in  the  areas  where  these  ^ 

i 

peculiarities  occur.  Such  situations  have  been  left  in  the  data  so  that  the  critical  j 

scientist  may  determine  how  much  faith  he  wishes  to  place  In  regional  means. 

(B)  Standard  deviations,  particularly  at  the  surface,  are  high  near  40°N,  i 

i 

175°W.  Notice  that  the  standard  deviations  drop  as  one  moves  eastward  along  latitude 

9 

40  and/or  with  Increasing  depth.  This  region  Is  In  the  vicinity  of  the  subarctic  front 
and  one  would  expect  high  variability  here.  Also,  since  it  occurs  In  the  westernmost  „ 
edge  of  the  study  region  It  may  be  ascribed  to  the  fact  that  this  region  Is  near  the  j 

extension  of  the  Kuroshio  and  Oyashio  confluence  and  thus  Is  an  area  where  high 
variability  Is  to  be  expected.  Visual  Inspection  of  the  data  that  have  gone  Into  the 
calculations  for  this  region  show  them  to  be  reasonable  with  respect  to  all  of  the 
verification  and  screening  of  standards  mentioned  above.  However,  the  data  appear 
to  characterize  two  rather  distinct  oceanic  states,  I.e.,  they  come  from  a subtropical 
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water  mass  or  a subarctic  water  mass.  Thus  the  estimates  of  the  mean  in  this  region 
may  be  somewhat  misleading  for  they  describe  a situation  which  physically  does  not 
exist.  Fortunately  this  effect  appears  confined  to  only  one  or  two  of  the  quadrangles 
in  the  vicinity  of  the  region  mentioned  and  does  not  affect  the  rest  of  the  mean 
field  calculations.  However,  the  user  of  these  data  is  forewarned  that  the  mean 
field  estimates  given  here  may  not  be  representative  of  the  actual  physical  conditions  : 
in  this  region. 

(C)  As  indicated  above,  a significant  portion  of  the  data  come  from  observa- 
tions at  the  ocean  weather  ships.  These  regions  of  high  observation  density  are 
clearly  apparent  in  the  Tables  shown  below.  It  is  important  to  note  that  the  estimates" 
of  the  mean  field  in  the  region  surrounding  the  weather  ships  have  in  no  way  been 
influenced  by  the  weather  ship  values.  On  the  other  hand,  if  one  wanted  to  construct 

an  areal  average  in  the  vicinity  of  the  weather  ships,  the  estimates  of  the  mean  i 

field  coming  from  those  boxes  with  a large  number  of  observations  probably  should  be  1 
given  more  weight  than  their  surrounding  neighbors. 

(D)  The  data  below  150  m come  mainly  from  XBTs.  It  is  clear  from  Figure  1 
that  these  data  are  abundant  only  after  1968.  Hence,  the  mean  field  at  greater 

depth  is  estimated  from  a relative  short  record.  This  depth  dependent  base  period  ^ 
for  construction  of  the  means  may  cause  problems  in  some  applications.  However,  it 
is  the  best  that  can  be  done  with  available  data.  This  problem  does  not  generally 
exist  in  the  vicinity  of  the  weather  ships  where  numerous  MBT  and  hydrographic  data 
extend  back  into  the  1950's. 

(E)  The  data  become  increasingly  sparse  with  greater  depth.  At  the  greatest  ‘ 
depths  the  data  are  presented  only  to  give  an  idea  of  the  total  coverage  to  date  rather 

j 

than  as  a valid  estimate  of  the  mean  field.  It  has  been  left  to  the  individual  ! 

scientist  to  judge  the  region/depth  below  which  the  estimates  become  meaningless. 
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Figure  1 Temporal  distribution  of  the  subsurface  temperature  data  by  month 
and  data  type  over  the  study  area. 

Figure  2 Seasonal  cycle  obtained  from  monthly  mean  data  at  depths  0 and  90  m 
for  two  typical  latitudes. 

Figure  3 Scatter  diagram  showing  the  relationship  between  the  sea-surface 

temperature  climatology  of  Namlas  (from  ship  Injection  temperatures), 
and  the  sea-surface  temperature  climatology  developed  In  this  work. 
Data  shown  for  all  months  In  the  study  region  defined  1r.  the  text. 


Figure  4 Distribution  of  the  number  of  occurrences  of  temperature  difference 
obtained  by  subtracting  the  BT  sea- surface  temperature  climatology 
from  the  sea- surface  temperature  climatology  of  Namlas.  The  bias  Is 
expected  from  the  results  of  earlier  workers. 
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- 1.12  Monthly  mean  temperatures  for  the  central  Pacific  at  the  sea  surface. 

- 2.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  30  m. 

- 3.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  60  m. 

- 4.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  90  m. 

- 5.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  120  m. 

- 6.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  150  m. 

- 7.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  200  m. 

- 8.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  250  m. 

- 9.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  300  m. 

- 10.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  400  m. 

- 11.12  Monthly  mean  temperatures  for  the  central  Pacific  at  a depth  of  500  m. 


THERMAL  STRUCTURE  OBS. 


Figure  2 

Seasonal  cycle  obtained  from  monthly  mean  data  at  depths  0 and  90  m for  two  typical  latitudes. 


Figure  3 

Scatter  diagram  showinq  the  relationship  between  the  sea  surface  temperature 
climatology  of  Namlas  (from  ship  Injection  temperatures),  and  the  sea  surface 
temperature  climatology  developed  In  this  work.  Data  shown  for  all  months  In 
the  study  region  defined  In  the  text. 
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